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SUMMARY 


A  procedure  is  described  for  obtaining  minimum-mass  designs  of  insulated 
composite  structural  panels.  The  panels  are  characterized  as  consisting  of  a 
structural  layer  and  an  insulation  layer.  The  panels  are  loaded  by  a  general 
set  of  inplane  forces  applied  to  the  structural  layer  and  a  time-dependent  tem¬ 
perature  applied  to  the  outer  surface  of  the  insulation  layer.  Temperature  and 
stress  histories  in  the  panel  are  given  by  closed-form  solutions,  and  optimiza¬ 
tion  of  the  insulation  and  structural  thicknesses  is  performed  by  nonlinear 
mathematical  programming  techniques.  Design  constraints  are  enforced  at  a 
finite  number  of  discrete  times  over  the  time  period  of  interest,  and  satisfac¬ 
tory  results  are  obtained  with  a  small  number  of  times. 

The  computerized  procedure  is  intended  for  preliminary  design  calculations 
to  evaluate  materials  for  specified  applications  and  to  perform  parameter  stud¬ 
ies.  A  set  of  design  calculations  is  presented  to  evaluate  the  efficiency  of 
the  following  eight  structural  materials  under  combined  heating  and  mechanical 
loads:  graphite/polyimide  (Gr/PI),  graphite/epoxy  (Gr/E) ,  boron/aluminum 

(B/Al)  ,  titanium  (Ti)  ,  aluminum  (Al)  ,  Rene*^  41  ,  carbon/carbon,  and  Lockalloy. 

A  study  is  also  performed  to  assess  the  effect  on  design  mass  of  the  intensity 
and  duration  of  heating  for  the  eight  materials.  Examinations  of  final  designs 
indicate  that,  for  sufficiently  high  mechanical  loads,  an  optimum  structure  has 
a  temperature  response  well  below  the  recommended  allowable  temperature  for  the 
material.  A  canparison  of  these  designs  with  those  assumed  to  operate  at  the 
allowable  temperatures  indicates  that  significant  mass  savings  may  be  attainable 
by  the  lower  temperature  operation. 


INTRODUCTION 

A  current  technical  challenge  in  the  design  of  high-speed  flight  vehicles 
is  the  development  of  optimum  structural  design  concepts  for  combined  heating 
and  mechanical  loading  applications.  A  concept  of  interest  for  applications  to 
atmospheric  reentry  vehicles  is  the  insulated  configuration  in  which  the  struc¬ 
ture  is  protected  from  heating  by  insulation  placed  between  the  airstream  and 
the  structure.  The  potential  benefits  of  structural  concepts  composed  of  fiber- 
reinforced  composite  materials  are  currently  being  evaluated  (ref.  1).  At  an 
early  stage  in  the  design  of  such  thermal-structural  concepts,  it  is  beneficial 
to  conduct  selected  design  studies  to  (1)  evaluate  candidate  structural  and 
insulation  materials  to  narrow  the  choice  of  materials  to  a  small  number  of  con¬ 
tenders  and  (2)  obtain  a  preliminary  estimate  of  the  optimum  combination  of 
structural  and  insulation  material.  An  emerging  philosophy  for  conducting  such 
studies  is  integrated  thermal-structural  analysis  and  design  (refs.  2  to  4)  in 
which  all  important  interactions  between  the  thermal  and  structural  responses 
of  the  structure  are  properly  accounted  for  at  every  stage  in  the  design  pro¬ 
cess.  Presented  herein  is  a  design  technique  to  define  the  best  combination 
of  insulation  and  structural  thicknesses  which  minimizes  the  total  mass  under 
a  specified  set  of  inplane  mechanical  loads  and  a  transient  heating  pulse. 


Previous  investigators  have  addressed  the  minimum-mass  design  of  insulated 
structures  (refs.  5  to  8) ,  but  consideration  has  been  limited  to  metallic  struc¬ 
tures,  uniaxial  loading,  and  a  restricted  class  of  heating  pulses.  These  sim¬ 
plifications  were  needed  since  the  classical  minimization  techniques  required 
extensive  manipulations  of  algebraic  equations. 

In  the  present  work  the  structure  considered  is  either  a  metal  layer  or  a 
ccxnposite  laminate  under  a  general  set  of  inplane  mechanical  loads  and  heating 
pulse.  This  generalization  is  facilitated  by  the  use  of  nonlinear  mathematical 
programming  techniques  to  generate  the  designs  systematically.  An  important 
aspect  of  the  present  approach  is  the  treatment  of  the  requirement  that  time- 
dependent  temperatures  and  stresses  be  less  than  specified  allowable  values  for 
all  appropriate  values  of  time.  This  aspect  has  received  relatively  little 
attention.  In  reference  9,  which  deals  with  the  optimization  of  an  ablating 
heat  shield,  constraints  on  time-dependent  temperatures  and  stresses  are 
replaced  by  time-integrated  averages  of  the  original  constraints.  The  inte¬ 
grals  are  carefully  formulated  to  avoid  constraint  violations,  but,  as  pointed 
out  in  reference  9,  the  integral  representation  of  the  constraints  leads  to  a^ 
smoothing  effect  which  causes  temperatures  and  stresses  to  be  somewhat  insensi¬ 
tive  to  changes  in  the  structural  size  parameters.  The  approach  used  in  the 
present  work  is  a  simple  time  slicing  in  which  the  constraints  are  enforced  at 
a  sufficiently  large  number  of  discrete  times  to  ensure  satisfaction  of  all  con¬ 
straints  throughout  the  time  period  of  interest.  (Since  the  completion  of  this 
work,  progress  has  been  reported  in  ref.  10  on  methods  to  track  one  or  more  most 
critical  tiroes  and  enforce  constraints  only  at  those  times.)  Analytical  solu¬ 
tions  are  obtained  for  the  temperature  history  and  stresses  in  the  structure. 

The  mass  minimizations  are  carried  out  by  use  of  the  general-purpose  optimizer 
AESOP  (ref.  11). 

The  analytical  design  procedure  is  applied  to  two  studies:  (1)  A  compara¬ 
tive  structural  efficiency  study  of  several  structural  materials  suitable  for 
high- temperature  structural  applications  -  graphite/polyimide,  graphite/ epoxy, 
boron/aluminum,  titanium,  Ren/  41 ,  aluminum,  carbon/carbon,  and  Lockalloy; 
and  (2)  a  study  of  the  effect  of  intensity  and  duration  of  heating  on  the 
required  mass  for  these  materials.  Effects  of  certain  simplifying  assumptions 
made  in  the  analysis  are  also  assessed  for  selected  final  designs.  The  purposes 
of  this  paper  are  to  describe  the  analytical  design  procedure  and  to  present 
the  results  of  the  studies  carried  out  with  the  procedure. 


SOBOLS 

The  values  are  given  in  both  SI  and  U.S.  Customary  Units.  The  calcula- 
tions  were  made  in  U.S.  Customary  Units. 

c  heat  capacity,  J/kg-°C  (Btu/lbrn-Op) 

ci,C2  heat  capacity  of  insulation  and  structural  material,  respectively, 
J/kg-°C  (Btu/lbm-OF) 

j),  ,Eo  Young's  modulus  in  fiber  and  transverse  directions,  respectively.  Pa 
(lbf/in2) 
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F  Tsai-Wu  failure  criterion 

F"!  =  1 /Xiji  +  1 /Xq 

F2  =  VYt  +  1/Yc 

Fll  =  -VXtXc 

F22  =  -VYtYc 

Fee  =  i/s2 

G  shear  modulus,  Pa  (Ibf/in^) 

g  constraint 

k  thermal  conductivity,  W/m~^C  (Btu/in-sec-^F) 

k*!  ,k2  thermal  conductivity  of  insulation  and  structural  material,  respec¬ 
tively,  W/m-^^C  (Btu/in-sec-^F) 

NLAYER  number  of  layers  in  laminate 

NTIMES  number  of  discrete  times 

Nj^,Ny,Nxy  stress  resultants,  kN/m  (Ibf/in.) 

thermal  forces,  kN/m  (Ibf/in.) 

RT  room  temperature 

S  shear  strength.  Pa  (Ibf/in^) 

T  temperature,  K  (^F) 

allowable  temperature,  K  (^F) 

Teq  outer  surface  temperature,  K  (^F) 

'^ego^'^peak^’^egf '1  parameters  of  temperature  pulse  (fig.  2),  K  (^F) 

^peak^^f  J 

T*  optimum  operating  temperature,  K  (^F)  (fig.  3) 

tO/t45,t9o  thickness  of  0^,  ±45°,  and  90°  plies,  respectively,  cm  (in.) 

t-|,t2  insulation  and  structure  thicknesses,  respectively,  cm  (in.) 

V  Von  Mises*  stress  criterion.  Pa  (Ibf/in^) 

W  mass,  kg  (Ibm) 
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Xt,Xc 


tensile  and  compressive  strengths  in  fiber  direction,  respectively. 
Pa  (Ibf/in^) 


Yt,Yc 


x,Y 


a-|,a2 


3 
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V-|  ,V2 


P 

Pi  rP2 


^1  r'^2'^12 
T 

T 


tensile  and  compressive  strengths  in  transverse  direction,  respec 
tively.  Pa  (Ibf/in^) 

laminate  coordinate  directions 

QQQj^dinates  through  depth  of  insulation  and  structure,  respectively, 
cm  (in.) 

coefficients  of  thermal  expansion  in  fiber  and  transverse  directions, 
respectively,  (°F"^) 

edge-fixity  coefficient 

ply  angle,  deg 

Poisson's  ratio  for  fiber  and  transverse  directions,  respectively 
mass  density,  kg/m^  (Ibm/in^) 

mass  density  of  insulation  and  structural  materials,  respectively, 
kg/m^  (Ibm/in^) 

allowable  stress.  Pa  (Ibf/in^) 

stress  components  in  laminar  coordinate  system.  Pa  (Ibf/in  ) 
time,  sec 

period  of  time  that  constraints  are  imposed;  monitoring  time,  sec 


CONFIGURATION  AND  LOADS 

The  configuration,  as  shown  in  figure  1 ,  consists  of  a  layer  of  insulation 
attached  to  a  structural  layer.  This  configuration  is  assumed  to  be  a  localized 
region  of  a  larger  structure  and  is  sufficiently  small  so  that  neither  the 
mechanical  loads  nor  the  temperature  vary  significantly  over  the  region.  Also 
the  region  is  assumed  to  be  sufficiently  far  removed  from  restrained  boundaries 
so  that  edge  effects  are  negligible.  The  structural  layer  is  subjected  to  a 
general  set  of  inplane  mechanical  loads  given  by  the  stress  resultants  N^/  Ny, 
and  Nxy,  which  are  assumed  to  be  time  independent.  The  structural  layer  is 
either  a  filamentary  composite  laminate  or  a  metal  layer.  No  mechanical  loads 
are  applied  to  the  insulation  layer.  The  outer  surface  of  the  insulation  is 
subjected  to  time^-dependent  heating,  which  is  characterized  by  a  temperature 

(T) .  In  the  present  work,  the  particular  form  of  is  the  curve  shown 

in^figure  2.  The  shape  of  the  curve  (generally  a  third-order  function  of  time 
for  T  ^  T£  and  constant  thereafter)  is  determined  by  the  quantities  '^eqo^ 

Tpeak'  ^peak'  ^f'  "^eqf* 
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ANALYTICAL  DESIGN  PROCEDURE 


Thermal  and  Structural  Analysis 

Analysis  of  the  temperature  history  is  based  on  one-dimensional  heat  trans¬ 
fer  through  the  depth  of  the  insulation-structure  configuration.  The  outer  sur¬ 
face  of  the  insulation  has  a  prescribed  temperature,  and  the  back  face  of  the 
structural  layer  is  assumed  to  be  an  adiabatic  surface.  The  temperature  history 
throughout  the  insulation  and  structure  is  obtained  by  an  analytical  solution 
described  in  appendix  A.  The  assumption  that  the  back  face  of  the  structural 
layer  is  an  adiabatic  surface  is  known  to  be  conservative.  Quantification  of 
the  effects  of  the  adiabatic  assumption  is  described  in  appendix  B,  where  this 
assumption  is  relaxed  by  permitting  radiation  from  the  back  face  of  the  struc¬ 
tural  layer,  and  temperatures  for  radiative  and  adiabatic  back-face  conditions 
are  compared. 

The  structural  analysis  is  based  on  calculating  stresses  and  strains  at  a 
point  in  the  structure  resulting  from  a  set  of  applied  forces  N^r  Ny,  and  N^y 
and  the  temperature  at  that  point.  This  approach  neglects  edge-boundary  effects, 
especially  those  due  to  thermal  stresses  arising  from  restraint  of  thermal 
expansion;  however,  thermal  stresses  due  to  interlayer  thermal  property  mismatch 
are  included.  The  applied  forces  are  assumed  to  be  constant  (i.e.,  time  inde¬ 
pendent)  .  The  structural  analysis  based  on  these  assumptions  is  described  in 
appendix  A.  Effects  of  neglecting  restrained  thermal  expansion  due  to  edge 
restraint  and  effects  of  time-dependent  loading  are  evaluated  in  appendix  B. 

In  the  analysis,  linear  representations  of  the  temperature  dependence  of  the 
mechanical  properties  are  used.  In  principle,  any  known  variation  of  mechanical 
property  with  temperature  may  be  input. 


Design  Objective  and  Constraints 

The  objective  of  the  present  procedure  is  to  determine  the  insulation 
thickness  ti  and  structural  thickness  t2  for  the  configuration  in  figure  1 
such  that  the  total  mass  is  a  minimum  and  the  structural  temperature  and 
stresses  do  not  exceed  prescribed  allowable  values  throughout  a  span  of  time 
T.  The  mass  is  given  by 


W  =  p-|t*|  +  p2t2 


(1) 


where  for  a  composite  laminate  with  0^,  ±45^,  and  90°  plies 


t2  =  to  +  t45  +  tgo 


(2) 


and  for  a  metal  structure  t2  is  the  metal  thickness.  The  temperature  con¬ 
straint  is 


5 


where  T  is  the  temperature  of  the  structural  layer  and 
temperature  of  the  structural  layer.  The  value  assigned  to 


is  the  allowable 
'  is  somewhat 


arbitrary  and  represents  the  maximum  service  temperature  of  the  material  for 
structural  applications.  In  the  case  of  filamentary  composites,  values  of 
are  generally  dictated  by  the  integrity  of  the  resin  material.  In  this  paper, 
values  of  T;,  are  based  on  current  usage  of  the  various  materials  in  high- 
temperature  applications  and  recommended  temperature  ranges  from  the  appropriate 
literature  (refs.  12  to  15). 

The  stress  constraint  for  each  layer  of  a  filamentary  composite  laminate 
is  expressed  by  a  modified  form  of  the  Tsai-Wu  failure  criterion  (ref.  16) 


F-jai  +  F2<^2  ^IT^I  *  ^22^2  ^66*^12 


(0  <  T  S  T) 


The  coefficients  F, ,  Fj,  F,,,  F22.  «ha  Fee 

of  the  allowable  stresses  for  the  composite  system  and  are  all  temperature 
dependent.  The  stress  constraint  for  a  metal  layer  is  given  by 


F  =  V/a-  ^  1 


(0  S  T  ^  t) 


where  the  Von  Mises'  stress  criterion  is 

/  2  2  ,  '  V''"'  (6) 

V  =  I  a-j  +  02  -  a-|a2  +  3a-|2l 

Before  describing  the  analytical  design  procedure,  it  is  instructive  to  describe 
the  trade-off  involved  in  designing  an  insulated  structure.  As  shown  in  fig¬ 
ure  3,  the  required  mass  of  the  load-carrying  structure  increases  with  operat- 
ing  temperature  of  the  structure,  primarily  because  of  degradation  of  mechanical 
properties.  The  required  mass  of  the  insulation  decreases  with  increasing  oper¬ 
ating  temperature  level.  The  total  mass  therefore  attains  a  minimum  for  some 
temperature  T*.  This  optimum  temperature  may  be  significantly  less  than 
allowable  temperature  indicated  by  T^.  It  is  shown  subsequently  that  the  dif¬ 
ference  between  T*  and  Ta  is  principally  a  function  of  the  mechanical  load. 


Design  Formulation 

The  design  problem  is  formulated  and  solved  as  a  nonlinear  mathematical 
programming  problem.  The  objective  function  is  given  by  equation  (1).  The 
design  variables  consist  of  ti  and  either  the  set  tQ,  ^45,  and  tgo  or 
to.  The  constraints  are  expressed  by  equation  (3)  and  either  equation  (  ) 
or  (5).  For  convenience  and  in  preparation  for  the  discretization  of  the  time 
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dependence  of  the  constraints,  equations  (3)  to  (5)  are  written  as  nondimen- 
sional  constraint  functions 

91 (T)  =  T(T)/Ta  -ISO  (0  S  T  S  T)  (7) 

gj+l(T)  =  Fj(T)  -  1  (0  S  T  S  T;  j  =  1,  2,  .  .  .  NLAYER)  (8) 

where  F  is  given  by  equation  (4)  or  (5) .  In  addition  to  the  behavior  con¬ 
straints  in  equations  (7)  and  (8),  minimum  gage  constraints  are  applied  to  the 
design  variables. 


Treatment  of  Constraints  on  Time-Dependent  Quantities 

In  developing  the  present  analytical  design  procedure,  special  attention 
was  focused  on  the  satisfaction  of  constraints  on  time-dependent  temperatures 
and  stresses  (eqs.  (7)  and  (8)).  Some  guidance  was  available  from  reference  9, 
which  treated  the  optimum  design  of  an  ablating  heat  shield.  In  reference  9 
the  constraints  on  time-dependent  quantities  were  replaced  by  time- integrated 
averages  of  the  original  constraints.  The  representations  were  carefully  for¬ 
mulated  to  exclude  violations  at  any  point  in  time  but,  as  pointed  out  in  ref¬ 
erence  9,  had  the  disadvantage  that  temperatures  and  stresses  were  somewhat 
insensitive  to  structural  size  parameters  because  of  the  smoothing  effect  of 
the  integral  representation  of  the  constraints.  Early  in  the  present  work, 
alternative  approaches  were  considered,  including  other  integral  representa¬ 
tions  (ref.  17)  and  ideas  for  tracking  one  or  more  most  critical  times  and 
enforcing  constraints  only  at  those  times.  The  approach  decided  upon  was  a 
simple  time  slicing,  whereby  constraints  were  enforced  at  a  sufficiently  large 
number  of  discrete  times  in  the  time  period  of  interest.  A  schematic  represen¬ 
tation  of  the  discrete-times  approach  for  temperature  constraints  is  shown  in 
figure  4. 

In  the  discrete-times  approach,  the  constraints  of  equations  (7)  and  (8) 
are  replaced  by 


91  =  T(Ti)/Ta  -ISO 


(i  =  1 ,  NTIMES)  (9) 


9NTIMES+i’ j 


Fj  (Ti)  -ISO 

(i  =  1 ,  NTIMES;  j  =  1 ,  2, 


.  NLAYER)  (10) 


Thus  the  total  number  of  constraints  is 


(1  +  NLAYER) (NTIMES) 
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Because  the  number  of  constraints  is  proportional  to  the  number  of  discrete 
times  NTIMES^  it  is  of  interest  to  determine  how  large  NTIMES  must  be  for 
typical  calculations.  The  answer  is  problem  dependent,  but  studies  carried  out 
and  documented  in  reference  4  and  in  a  subsequent  section  of  this  paper  indicate 
that,  for  the  types  of  insulated  structures  under  the  types  of  heating  histories 
considered  herein,  only  a  few  (three  or  four)  discrete  times  are  needed. 


Sizing  Technique 

The  computer  program  used  to  perform  the  optimizations  is  a  general-purpose 
optimizer  denoted  AESOP  (ref.  11),  which  accounts  for  constraints  using  an  exte¬ 
rior  penalty  function  approach  (ref.  17).  The  AESOP  code  contains  several  opti¬ 
mization  algorithms  that  can  be  selected  in  various  combinations  by  the  user. 

The  approach  used  here  was  a  combination  of  adaptive  search  and  pattern  search. 
These  algorithms  are  described  in  reference  11. 

APPLICATIONS  AND  DISCUSSION 
Overview 

In  this  section  the  design  procedure  is  used  to  carry  out  design  studies. 
First  a  preliminary  study  is  performed  to  assess  the  minimum  number  of  discrete 
times  needed  to  obtain  minimum-mass  designs.  Next  the  design  procedure  is 
applied  to  an  efficiency  study  for  eight  structural  materials  over  a  wide  range 
of  loads  and  two  heating  conditions  representative  of  high-speed  flight  of 
winged  reentry  vehicles  (refs.  18  and  19) .  The  materials  in  the  study  include 
graphite/epoxy,  graphite/polyimide,  boron/aluminum,  aluminum,  titanium, 

Rene  41,  carbon/carbon,  and  Lockalloy.  A  study  is  then  carried  out  to  assess 
the  effect  on  minimum  mass  of  heating  duration  and  intensity  for  the  eight 
materials.  Finally  results  are  presented  to  assess  the  effects  of  certain  sim¬ 
plifying  assumptions  made  in  the  analysis. 


Materials  Used  In  Studies 

The  eight  structural  materials  selected  for  the  studies  are  listed  in 
table  I  along  with  averaged  values  used  for  thermal  conductivity,  specific 
heat  capacity,  and  density,  as  well  as  the  allowable  temperatures  used  in  the 
studies.  Properties  are  also  given  for  the  insulation  which  is  LI-900  reus¬ 
able  surface  insulation.  Mechanical  properties  for  the  structural  materials 
were  obtained  from  references  4  and  12  to  15  and  are  listed  in  table  II.  Min¬ 
imum  gages  considered  are  given  in  table  III . 


Discrete-Times  Convergence  Study 

The  purpose  of  this  section  is  to  show  approximately  how  many  discrete 
times  are  needed  for  a  typical  design  problem. 
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The  study  is  carried  out  for  a  Gr/E  panel  with  an  outer  surface  temperature 
history  (fig.  2)  characterized  by  Tpeak  =  K  (1500O  F)  ,  Tp^ai^  =  700  sec, 

Tf  =  1500  sec,  and  Tgqf  =  297  K  (75^  F)  and  loads  of  =  Ny  =  -700  kN/m^ 
(-4000  Ibf/in.)  and  =  525  kN/m  (3000  Ibf/in,) .  The  monitoring  time  f 

was  6000  sec.  Designs  were  obtained  for  3  to  1 5  time  slices.  Final  mass  is 
plotted  against  the  number  of  slices  in  figure  5.  The  plot  indicates  that  only 
three  or  four  discrete  times  are  necessary  to  give  satisfactory  results  in  the 
sense  that  increasing  the  number  of  discrete  times  does  not  improve  the  result. 
Figure  6(a)  shows  the  temperature  and  stress  response  (in  the  45°  ply)  in  the 
final  design.  The  response  of  temperature  is  most  important  to  this  design 
(i,e,,  only  the  temperature  constraint  is  critical  for  T  >  0) ,  Observe  that 
the  temperature  response  is  a  smooth  curve  with  a  single  local  maximum.  Thus 
only  a  small  number  of  equally  spaced  time  slices  are  needed  to  define  properly 
the  critical  temperature.  If  a  more  severe  or  complicated  heating  pulse  is 
applied,  the  temperature  response  would  still  be  smooth  and  gradual.  This 
assertion  is  based  on  previous  analyses  (refs,  4  and  20)  and  occurs  because 
the  insulation  has  a  smoothing  (dampinglike)  effect  on  the  thermal  response  of 
the  structure. 

The  stress-response  curve  is  somewhat  more  complicated  than  the  temperature 
response  but  again  only  has  a  single  local  maximum.  Since  stress  is  not  criti¬ 
cal  for  T  >  0,  it  is  not  clear  from  this  example  what  effect  the  stress 
response  has  on  the  number  of  discrete  times  required.  To  investigate  this 
effect,  consider  the  Gr/E  panel  with  higher  loads  of  =  Ny  =  -1.05  MN/m 
(-6000  Ibf/in.)  and  Nxy  =  700  kN/m  (4000  Ibf/in.),  in  which  stresses  are  crit¬ 
ical  for  T  >  0  (fig,  6(b)).  The  plot  of  final  mass  against  the  number  of  time 
slices  is  shown  in  figure  7  and  indicates  that,  despite  the  criticality  of 
stress,  only  three  or  four  time  slices  are  needed  for  a  satisfactory  answer. 


Thermal-Structural  Efficiency  Study 

This  study  is  a  comparative  evaluation  of  the  performance  of  the  eight 
materials  over  a  wide  spectrum  of  load  levels  for  two  heating  conditions. 

Except  for  Rene  41  and  carbon/carbon,  each  material  is  protected  by  a  layer 
of  insulation.  The  first  (low)  heating  condition  is  characterized  by 
Tpeak  “  '^089  K  (1  500°  F)  ,  and  the  second  (high)  heating  condition  has  a  value 
of  Tpgak  ^  K  (2300°  F)  .  In  both  cases  =  700  sec  and 

Tf  =  2000  sec.  In  all  calculations  Ny  =  and  N^y  =  (2/3) N^.  The  results 
for  the  low-heating  condition  are  presented  in  figure  8.  At  high-load  levels 
Gr/PI  gives  the  lowest  mass,  followed  by  B/Al.  Neither  Rene  41  nor  carbon/ 
carbon  are  competitive  at  high-load  levels.  For  Njj  between  175  and  525  kN/m 
(1000  and  3000  Ibf/in.)  Rene  41  (uninsulated)  and  Lockalloy  give  the  lowest 
mass,  and  for  N^  less  than  175  kN/m  (1000  Ibf/in.)  carbon/carbon  gives  the 
lowest  mass.  These  results  suggest  possible  benefits  of  these  unprotected 
materials  for  lightly  loaded  regions,  such  as  control  surfaces.  The  curves  for 
Lockalloy  and  aluminum  have  a  corner  or  discontinuity  that  represents  a  change 
in  the  character  of  the  final  design  -  designs  to  the  left  of  the  corner  are 
operating  at  the  allowable  temperature  and  allowable  stress,  whereas  designs 
to  the  right  of  the  corner  are  operating  at  the  allowable  stress  but  with  tem¬ 
peratures  less  than  the  allowable  value. 
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For  the  high-heating  case  in  which  Tp^a]^  =  1  533  K  (2300^  F)  (fig-  9), 
Gr/PI  again  gives  the  lowest  mass  for  the  highest  load  values  (above  about 
1.05  MN/m  (6000  Ibf/in.)).  For  loads  between  700  kN/m  and  1.05  MN/m  (4000  and 
6000  Ibf/in.),  B/Al  designs  have  the  lowest  mass.  Lockalloy  gives  the  lowest 
mass  designs  for  between  262  and  700  kN/m  (1500  and  4000  Ibf/in.),  and 

carbon/carbon  gives  the  lowest  mass  for  loads  less  than  262  kN/m  (1500  Ibf/in.). 
Based  on  these  results  uninsulated  car bon/ car bon  and  insulated  Lockalloy  again 
seem  to  have  significant  potential  for  application  to  lightly  loaded  regions, 
and  for  more  highly  loaded  regions  such  as  the  wing,  insulated  Gr/PI,  B/Al, 

Gr/E,  and  titanium  seem  to  be  worthy  of  consideration.  These  results  are  based 
on  a  simplified  model  and  consider  only  strength  and  temperature  requirements. 


Characteristics  of  Final  Designs 

The  main  characteristics  of  selected  designs  in  figures  8  and  9  are  sum¬ 
marized  in  tables  IV(a)  to  IV(1)  .  Tables  IV(a)  to  IV(f)  have  results  for  the 
lower  heating  case  where  Tp^ak  =1089  K  (1500®  f)  for  values  of  N^  from 
0  to  -2.63  MN/m  (-15  000  Ibf/in.)  in  increments  of  525  kN/m  (3000  Ibf/in.). 
Tables  IV (g)  to  IV(1)  have  corresponding  results  for  the  higher  heating  case 
where  Tp^ak  =  1 ^  (2300^  F) .  The  tables  contain  for  each  material  the  final 
mass,  maximum  temperature  response  in  the  final  design  and  the  time  it  occurs, 
and  maximum  stress  ratio  and  the  time  it  occurs.  The  tabulations  are  presented 
to  provide  check  results  for  other  investigators  and  to  show  some  features  of 
the  designs.  For  the  lowest  loads,  the  structures  operate  at  their  allowable 
temperatures  and  the  design  problon  is  one  of  providing  minimum-mass  designs 
for  temperature  requirements.  As  higher  loads  are  considered,  strength  con¬ 
straints  become  important  and  a  trade-off  mechanism  becomes  operative,  whereby 
high-temperature  operation  with  less  insulation,  lower  allowable  stress,  and 
a  heavier  structure  is  balanced  against  lower  temperature  operation  with  more 
insulation,  higher  allowable  stress,  and  a  lighter  structure.  For  each  insu¬ 
lated  material,  a  load  value  is  eventually  reached  where  the  trade-off  favors 
lower  temperature  operation  to  maintain  a  sufficiently  high  allowable  stress. 

As  suggested  by  figure  3,  the  structural  temperature  in  such  designs  is  below 
the  allowable  temperature  of  the  material  and  the  designs  are  governed  primarily 
by  strength  constraints. 

In  some  cases,  the  temperature  in  the  final  design  is  well  below  the  allow¬ 
able  temperature.  For  example,  Lockalloy  under  a  load  of  N^  =  -2.63  MN/m 
(-15  000  Ibf/in.)  and  Tp^ak  ^  (2300^  F)  attains  a  maximum  temperature 

which  is  only  49  percent  of  its  allowable  temperature.  This  effect  points  out 
that  contrary  to  conventional  design  practice,  the  minimum-mass  insulated  heat¬ 
sink  structure  is  not  necessarily  one  in  which  the  structure  operates  at  its 
allowable  temperature.  This  observation  was  previously  made  for  an  ablation 
panel  in  reference  9.  To  bring  out  the  significance  of  this  observation,  con¬ 
sider  a  set  of  Lockalloy  designs  for  T  =  1533  K  (2300^  F)  at  the  load  range 
of  0  to  -2.63  MN/m  (-15  000  Ibf/in.),  as  described  in  table  IV.  A  set  of  corre¬ 
sponding  designs,  based  on  the  conventional  practice  of  assuming  operation  at 
the  allowable  temperature,  is  generated  for  comparison.  The  latter  designs  are 
obtained  by  first  sizing  the  structural  thickness  from  the  applied  mechanical 
loads  and  the  allowable  stress  at  700  K  (800^  F) ,  then  fixing  the  structural 
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thickness,  and  then  carrying  out  design  calculations  to  size  the  insulation 
thickness  required  to  maintain  this  temperature.  Design  masses  from  this  con¬ 
ventional  procedure  and  the  present  method  are  compared  in  figure  10.  Since 
the  structure  has  a  peak  temperature  response  at  or  close  to  the  allowable  tem¬ 
perature,  there  is  little  difference  between  the  designs  at  loads  below  about 
525  kN/m  (3000  Ibf/in.) .  Significant  differences  occur  at  higher  loads.  For 
example,  at  2.10  MN/m  (12  000  Ibf/in.)  the  conventional  design  has  40  percent 
more  mass  than  the  present  design,  and  at  2.63  MN/m  (15  000  Ibf/in.)  the  conven¬ 
tional  design  has  47  percent  more  mass.  These  results  suggest  that  design  of 
structures  for  thermal  applications  should  be  based  on  simultaneous  considera¬ 
tion  of  temperature  and  structural  requirements. 


Effect  of  Heating  Duration  and  Intensity 

This  section  describes  a  study  of  the  effects  of  duration  and  intensity 
of  heating  on  minimum  mass.  Designs  are  obtained  corresponding  to  high  and  low 
values  of  intensity  and  duration  of  heating  associated  with  winged  reentry 
vehicles  (ref.  18),  as  described  by  the  parameters  Tpgak  ^^d  Xf,  respec¬ 
tively.  Values  used  for  Xf  were  1500  sec  and  2000  sec,  and  values  for 
Tpeak  1 089  K  (1500°  F)  and  1533  K  (2300°  F)  . 

Results  of  the  study  are  presented  in  table  V.  Both  intensity  and  duration 
of  heating  have  significant  effects  on  the  designs.  For  example.  Increasing 
Xf  for  B/Al  with  Tpgai^  =  1533  K  (2300°  F)  led  to  an  increase  in  required  mass 
of  18  percent.  Increasing  both  Tpeak  ^f  same  case  led  to  a  mass 

increase  of  41  percent. 


Effects  of  Simplifying  Assumptions 

This  section  summarizes  the  effects  of  certain  assumptions  made  in  the 
analysis;  namely,  neglect  of  edge  effects  due  to  restraint  of  thermal  expansion, 
neglect  of  time  dependence  of  the  mechanical  loads,  and  neglect  of  heat  loss 
from  the  back  face  of  the  structural  layer.  Details  are  given  in  appendix  B. 

Restrained  thermal  expansion.-  Calculations  are  carried  out  for  Gr/E, 
titanium,  and  Lockalloy.  Figure  11  shows  a  plot  of  design  mass  as  a  function 
of  an  edge-fixity  coefficient  3  which  varies  between  0  and  1 .  A  value  of 
3=1  corresponds  to  completely  restrained  expansion,  and  3=0  corresponds 
to  unrestrained  expansion.  Results  indicate  that  the  effects  of  restrained 
thermal  expansion  are  negligibly  small  for  Gr/E,  are  moderate  for  titanium  (a 
22-percent  increase  in  design  mass  as  3  is  increased  from  0  to  1 ) ,  and  large 
for  Lockalloy  (an  83-percent  increase  as  3  is  increased  from  0  to  1 ) . 

Time-dependent  loading.-  To  assess  the  effects  of  time  dependence  of  the 
mechanical  loads,  design  calculations  are  repeated  for  two  selected  cases: 

(1)  Graphite/epoxy  with  a  maximum  load  =  -1.05  MN/m  (-6000  Ibf/in.);  and 

(2)  titanium  with  a  maximum  load  of  -1.58  MN/m  (-9000  Ibf/in.).  The  mechanical 
loads  varied  with  time  according  to  figure  12,  which  is  representative  of  load 
against  time  during  reentry  (ref.  21).  Designs  are  compared  with  those  based 
on  constant  loading.  Results  are  given  in  table  VI,  which  shows  that  time- 
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dependent  loading  has  a  very  small  effect  on  the  final  design  mass.  Observe 
however  that  in  the  designs  for  time-dependent  loads,  the  maximim  stress  ratio 
occurred  at  the  times  of  maximum  load  (1600  sec)  for  both  materials  rather  t  an 
times  of  peak  temperature  (2800  sec  for  Gr/E  and  2000  sec  for  titanium)  when 
loads  were  constant.  There  was  no  appreciable  effect  on  the  temperature 
responses  for  either  material  as  a  result  of  time-dependent  loads. 

Heat  transfer  from  back  face  of  structure.-  To  assess  the  effect  of 
neglecting  heat  transfer  from  the  back  face  of  the  structural  layer,  thermal 
analyses  of  selected  designs,  which  included  radiative  heat  transfer  from  the 
back  face,  were  carried  out.  A  f inite-elCTient  thermal  model  was  formulated  by 
using  the  SPAR  thermal  analyzer  (ref.  22) .  In  the  SPAR  model  the  back  face  of 
the  structural  layer  was  permitted  to  radiate  to  a  medium  at  room  temperature. 
Temperature  histories  were  calculated  for  the  eight  final  designs  in  table  IV(d) 
with  the  adiabatic  back  face  and  the  radiating  back  face.  Table  VII  contains 
a  comparison  of  the  peak  temperatures  for  the  two  conditions.  The  largest 
change  in  peak  temperature  is  33  percent  for  titanium.  This  significant  temper¬ 
ature  difference  does  not  necessarily  translate  into  a  large  effect  on  design 
mass.  AS  an  example,  the  increase  in  strength  of  titanium  from  548  K  to  451  K 
(527°  F  to  352°  F)  is  about  16  percent.  If  the  structural  mass  of  the  titanium 
design  is  assumed  to  be  decreased  by  16  percent  while  the  insulation  mass 
remains  unchanged,  the  total  required  mass  is  decreased  by  about  11  percent. 


CONCLUDING  REMARKS 

A  procedure  is  described  for  obtaining  minimum-mass  designs  of  insulated 
composite  or  metal  structural  panels  subjected  to  transient  heating  and  mechani¬ 
cal  loads.  The  procedure  uses  nonlinear  mathematical  programming  techniques, 
and  analytical  solutions  are  used  for  temperature  histories  and  stresses  in  the 
structure.  An  important  aspect  of  the  present  problem  is  the  need  to  enforce 
constraints  on  transient  temperatures  and  stresses.  In  the  present  procedure, 
the  constraints  are  enforced  at  a  finite  number  of  time  slices  or  discrete  times 
over  the  time  period  of  interest,  and  satisfactory  results  are  obtained  with  a 
small  number  of  times . 

To  simplify  the  analysis,  edge  effects  due  to  restrained  thermal  expansion 
are  neglected,  constant  (time  independent)  mechanical  loads  are  assumed,  and 
heat  loss  from  the  back  face  of  the  structure  is  neglected.  The  errors  in 
design  mass  introduced  by  these  assumptions  are  examined.  Neglecting  edge 
effects  can  lead  to  sizable  errors;  neglecting  heat  loss  from  the  back  face  has 
a  smaller  but  significant  effect;  and  time  dependence  of  loading  has  a  negligi¬ 
ble  effect  on  the  designs. 

The  design  procedure  is  used  to  perform  structural  efficiency  studies  for 
eight  materials  for  two  representative  values  of  peak  external  heating  and  a 
range  of  mechanical  load  levels.  Materials  considered  are  graphite/polyimide, 
graphite/epoxy,  boron/aluminum,  titanium,  aluminum,  Rene  41 ,  carbon/carbon, 
and  Lockalloy.  All  except  Rene  41  and  carbon/carbon  have  a  layer  of  reusable 
surface  insulation.  For  the  lower  heating  condition  graphite/polyimide  gave 
the  lowest  mass  except  at  the  lowest  load  levels,  whereas  Rene  41  and  carbon/ 
carbon  gave  the  lowest  mass.  For  the  higher  heating  condition  graphite/ 
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polyimide  again  gave  the  lowest  mass  for  high  and  moderate  loads,  whereas 
Lockalloy  and  carbon/carbon  gave  the  lowest  mass  at  low-load  levels.  Addi¬ 
tional  calculations  indicate  that  intensity  and  duration  of  heating  can  have 
a  substantial  effect  on  design  mass. 

Examinations  of  the  characteristics  of  the  final  designs  indicate  that  for 
sufficiently  high  mechanical  loads,  the  optimum  structure  has  a  temperature 
response  substantially  less  than  the  recommended  allowable  temperature  of  the 
material.  This  result  is  contrary  to  the  usual  engineering  practice  of  basing 
a  design  on  the  higher  temperature  operation.  Conparisons  of  several  Lockalloy 
designs  obtained  by  the  present  procedure  with  those  based  on  operation  at 
higher  temperatures  indicate  that  significant  mass  savings  may  be  attainable  by 
lower  temperature  operation. 


Langley  Research  Center 

National  Aeronautics  and  Space  Administration 
Hampton,  VA  23665 
October  11,  1979 
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APPENDIX  A 


THERMAL  AND  STRUCTURAL  ANALYSIS 


TEMPERATURE  HISTORY 

The  solution  for  the  temperature  history  is  obtained  by  using  the  analysis 
of  heat  transfer  in  a  two-layered  configuration  given  in  reference  23.  The  pur 
pose  of  this  appendix  is  to  outline  the  approach  used  in  reference  23  and  spe¬ 
cialize  that  general  result  for  the  present  problen  (fig.  1).  The  governing 
equations,  boundary  conditions,  and  initial  conditions  are 


9T2 

^2  Z2.0 


0 


T]  (Z-,  ,0)  =  T2(Z2,0)  =  To 


(A6) 


(A?) 


Equations  (Al)  and  (A2)  govern  the  temperatures  in  the  insulation  and  structural 
layers,  respectively.  Equations  (A3)  and  (A4)  define  the  continuity  of  heat 
flux  and  temperature  at  the  insulation-structure  interface.  Equation  (A5)  is 
the  specification  of  the  temperature  history  at  the  upper  surface  of  the  insula¬ 
tion.  Equation  (A6)  states  that  no  heat  is  lost  through  the  back  face  of  the 
structural  layer.  Equation  (A7)  defines  the  initial  conditions  on  temperatures. 
As  pointed  out  in  reference  8,  the  assumption  that  no  heat  is  transferred  from 
the  back  face  of  the  structural  layer  is  conservative,  especially  when  radiation 
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from  the  back  face  to  internal  heat  sinks  represents  a  significant  heat  loss. 

In  such  cases y  however,  it  would  probably  be  necessary  to  insulate  the  back  face 
of  the  structure  to  protect  internal  equipment  from  excessive  heating  by  thermal 
radiation.  Further  discussion  of  the  effects  of  this  assumption  is  presented 
in  appendix  B.  Solutions  to  the  equations  are  obtained  in  reference  23.  The 
expression  for  the  average  temperature  in  the  structural  layer  is  of  interest 
for  the  current  problem.  This  solution  is 


T(T)  .  Te,(0)  +  [To  -  Teq(0)J  Y 


*  r"  - 


1"  .y2|p.o, 


where  Yn  root  of 


35  tan  Y  tan  -  =  1 

3 


^■[^2P2^2 
^2^^ Pi 


(A10) 


(AH) 


Yn  /  1  \  Yn  Yn 

Hn  =  —  (3?  +  -  )  cos  Yn  sin  —  +  (1  +  ^)  sin  Yn  cos  — 

2  \  p/  p  p 


(A12) 


and  dimensionless  time  is 


C2P2^2 


(A13) 
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In  the  present  work  the  outer  surface  temperature  Teq(T)  is  represented  as 
(see  fig.  2) 


Teq(T)  =  aeq  +  bgqT  +  CgqX^  +  dgqX^ 


T  S  Tf)  (AT  4) 


“  "^eqf 


(T  >  Tf)  (A15) 


The  coefficients  in  equation  (A14)  are  expressed  in  terms  of  the  parameters  of 
the  curve  in  figure  2.  The  conditions  for  evaluating  the  coefficients  are  that 
the  curve  of  equation  (A14)  passes  through  the  three  points  noted  on  the  curve 
and  has  a  zero  slope  at  T  =  Tpeak*  The  results  are 


Teqf  “  Teqo  "^eqo 


^peak 


^f^Tpeak  Tf) 


Teqo  “  Tpeak 
Tpeak 


^Tpeak^eq 


2 

beq  =  ”2Xpga]^Ceq  -  3Tpgai5deq 


^eq  “  ^eqo 


The  coefficients  bgq,  Cgq,  and  dg^  are  nondimens ionali zed  by  defining 


where 


^eq  ”  ^eq/^ 
Yeq  ~  ^eq/^^ 
^eq  =  eq/R^ 


2 

C2P2^2 


Substituting  equations  (Al4)  and  (A15)  into  equation  (A8)  gives  the  following 
solutions: 
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Teq(O)  +  SeqP  +  YeqP^  [j^O  ^ 


sin  Yn  2 
- ^  e-^nP 


y  -  J"-)  .  *  Y^P  -  1 

4;  YnHn  ^  ^ 

n=l  'n  'n 


,  3.,,  ^  ^  Y  1(,  -  .-VnP 

y2  y'^  Y®^ 

'n  'n  'n 


T  ^  Tf)  (A16) 


T  =  Torrf  + 


uu 

[to  -  T^q(O)]  2 


!!1Z2  e-'pPf 


.  y  2!J[;  ^  jn^Pf-P)  .  J^P 


+  _[f2  e-YnP  .  Jn(P£-P)|',  . 

■v4 


-  36eq(^ 


!£!  e'^n(Pf-P) 


+  LLe^n(Pf-P) 


(T  >  Tf)  (AT  7) 


where 


^2 

Pf  =  - 2 

C2P2t2 


In  the  evaluation  of  equations  (AT  6)  and  (AT  7) ,  six  terms  of  the  series  have  pro¬ 
vided  sufficiently  accurate  temperatures.  In  the  general  case  where  both  t] 
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and  t2  are  nonzero,  the  roots  of  equation  {A9)  are  extracted  numerically  by 
a  systematic  search.  For  the  special  case  where  t]  =0,  that  is,  the  insula¬ 
tion  layer  is  emitted,  the  roots  of  equation  (A9)  are  given  by 


-  n  -  -J. 


{A18) 


and  for  this  case  equation  (A12)  reduces  to 


i  ’i 


(A19) 


STRESS  ANALYSIS 

The  stresses  in  the  laminate  are  computed  from  elementary  thermoelastic 
lamination  theory  (ref.  24) •  The  constitutive  equations  for  a  balanced  symmet¬ 
ric  laminate  under  inplane  mechanical  and  thermal  loading  are  written 


Ny  + 


An 

A12 

A]  2 

^22 

0 

0 

0  e 


0  /e 


-  yy 

^66  ®xy 


{A20) 


where  N^,  Ny,  and  Nxy  are  constant  (i.e.,  time  independent)  applied  mechani 
cal  forces  per  unit  width  of  the  laminate  and  An,  A]2f  A22r  and  A55  are 
laminate  stiffnesses,  given  by 


NLAYER 


(A21a) 


NLAYER 


(A21b) 
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NLAYER 


^22  ==  7  Q22H 


(A21c) 


NLAYER 


=  ^  Qeeti 


(A21d) 


where  ti  is  the  thickness  of  the  ith  layer.  The  thermal  forces  Ntx  and 
N^y  are  given  by 


NLAYER 


Ntx  =  ^  [(Qn^x  +  Ql2«y  +  Ql6«XY)t]i(T  -  Tj-ef) 


(A22) 


NLAYER 


Nty  =  ^  [jQl2°‘x  Q22“y  Q26“xy)  i  “  '•^ref) 


(A23) 


where  T^-gf  is  the  stress-free  temperature  and 


ax  =  +  n^a2 


a„  =  n'^a-i  +  m'‘a2 


ttxv  =  2mn(a-|  -  02) 


(A24) 


In  these  equations,  m  =  cos  9,  n  =  sin  9,  9  is  the  lamina  ply  angle,  and 


Qn  =  Qnm^  +  2(Q-|2  +2Q66)n^^  +  022"' 


(A25a) 


Q22  =  Qll"^  2(Qi2  +2066)"^™^  022™^ 


(A25b) 
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Q12  =  (Qn  +  Q22  -  4Q66)nV  +  Qi2(n'*  +  m^)  (A25c) 

Q16  =  (Qn  -  Q12  -2Q66)™n3  +  (Qi2  -  Q22  +  2Q66)n3m  {A25d) 

Q26  =  (Qn  "  Q12  -2Q66)'’^™  +  (Q12  -  Q22  +  2Q66)nm3  {A25e) 

Q66  =  (Qn  -  Q22  -  2Qi2  -  2Q66)in2n2  +  Qee  (n^  +  m^)  {A25f) 


The  stresses  in  equations  (A28)  neglect  edge  effects,  particularly  effects  of 
restrained  thermal  expansion  due  to  edge  supports*  Estimates  of  the  importance 
of  such  thermal  stresses  are  considered  in  appendix  B. 
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EFFECTS  OF  RESTRAINED  THERMAL  EXPANSION,  TIME- DEPENDENT  LOADS,  AND 
BACK-FACE  THERMAL  BOUNDARY  CONDITION 

The  purpose  of  this  appendix  is  to  quantify  effects  of  certain  simplifi¬ 
cations  employed  in  the  analysis.  These  simplifications  include  neglect  of  edge 
effects  due  to  restrained  thermal  expansion,  neglect  of  time  variations  of  the 
applied  mechanical  loads,  and  neglect  of  heat  loss  from  the  back  face  of  the 
structural  layer. 


RESTRAINED  THERMAL  EXPANSION 

The  stress  analysis  used  in  the  design  procedure  and  described  in  appen¬ 
dix  A  is  based  on  stresses  and  strains  at  a  point  and  thus  does  not  incorporate 
edge  effects,  particularly  effects  of  thermal  stresses  due  to  restraint  of  ther¬ 
mal  expansion  by  edge  fixity.  The  following  analysis  is  an  attempt  to  account, 
in  an  approximate  manner,  for  the  effects  on  design  mass  of  thermal  stresses 
due  to  edge  fixity.  If  the  temperature  of  a  panel  is  raised  while  thermal 
expansion  is  prevented,  the  panel  develops  forces  given  by 


Nx  =  -Ntx 

Ny  =  -Nty  > 
Nxy  =  0 


(B1) 


when  Nj-x  and  Nj-y  are  thermal  forces  defined  in  appendix  A.  For  the  special 
case  of  an  isotropic  layer 


Eoh 

^X  ~  Ny  =  -  (T  -  Tj.gf ) 

1  -  V 


(B2) 


where  E  is  Young's  modulus,  a  is  the  coefficient  of  thermal  expansion, 

V  is  Poisson's  ratio,  and  h  is  the  thickness  of  the  layer.  Equations  (Bl) 
correspond  to  complete  restraint  of  inplane  displacement.  For  completely  unre¬ 
strained  expansion,  the  forces  N^  and  Ny  are  zero.  This  circumstance  sug¬ 
gests  the  introduction  of  an  edge-fixity  coefficient  g,  in  which  thermal 
stresses  due  to  edge  restraint  are  given  by 


(B3) 


Nx  =  -3Ntx 

Ny  =  -gNty 
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where 


0  ^  6  ^  1 


Completely  fixed  edges  correspond  to  B  -  1 ,  and  free  edges  correspond  to 
R  =  0.  The  analysis  in  appendix  A  is  modified  to  incorporate  equations  (B3) 
byreplacing  Ntx  by  (1  -  3)Ntx  and  Nfcy  by  (1  -  3)Nty  inequation  (A20) . 

Evaluation  of  thermal-stress  effects  is  carried  out  for  three  materials: 
Lockalloy,  titanium,  and  graphite/epoxy.  These  materials  were  chosen  because 
their  properties  result  in  high,  intermediate,  and  low  values  of 
(See  app.  A  and  eq.  (B2) .)  The  results  are  shown  in  figure  11,  wherein  the  opti¬ 
mum  mass  for  the  three  materials  is  plotted  as  a  function  of  3.  Designs  corre¬ 
spond  to  Nx  =  -1.58  MN/m  (-9000  Ibf/in.),  Tpeak  =1089  K  (1500  F) ,  and 
Tf  =  2000  sec.  The  designs  for  3=0  are  listed  in  table  iy(d) .  Figure 
indicates  that  thermal  stress  due  to  edge  fixity  may  have  a  significant  effect 
on  mass.  This  effect  is  apparent  for  Lockalloy,  wherein  the  mass  corresponding 
to  3=1  is  83  percent  higher  than  the  mass  when  3=0.  The  change  for  tita¬ 
nium  is  22  percent  and  essentially  zero  for  graphite/epoxy.  The  importance  of 
the  thermal  stress  term  is  situational  dependent.  These  results  suggest  the 
need  to  extend  the  design  procedure  to  finite  panels,  wherein  edge  effects  would 
naturally  be  included.  Another  benefit  from  extending  the  current  procedure  to 
finite  panels  is  to  enable  buckling  requirements  to  be  properly  included  in  the 
design  of  insulated  panels. 


TIME-DEPENDENT  LOADING 

Design  calculations  presented  and  discussed  in  the  main  text  are  based  on 
constant  applied  mechanical  loads.  In  an  actual  application,  such  as  reentry, 
the  loads  vary  with  time.  This  section  describes  calculations  to  examine 
designs  of  time  dependence  of  the  loads.  Figure  12  depicts  a  representative 
plot  of  amplitude  of  force  (load  factor)  against  time  during  reentry  (ref .21 ) . 
AS  indicated  in  figure  12,  all  three  load  components  are  assumed  to  have  the 
same  time  variation. 

Calculations  are  carried  out  for  two  selected  designs:  Graphite/epoxy  with 
a  maximum  load  of  Nx  =  -1-05  MN/m  (-6000  Ibf/in.)  and  titanium  with  a  maximum 
load  of  -1.58  MN/m  (-9000  Ibf/in.).  Results  given  in  table  VI  show  that  time- 
dependent  loading  has  a  very  small  effect  on  the  final  design  mass.  I"  bhe 
designs  for  time-dependent  loads,  the  maximum  stress  ratio  occurred  at  the  times 
of  maximum  load  (1600  sec)  for  both  materials  rather  than  at  the  times  of  peak 
temperature  (2800  sec  for  Gr/E  and  2000  sec  for  titanium)  when  loads  were  con¬ 
stant.  There  was  no  appreciable  effect  on  the  temperature  responses  for  either 
material  as  a  result  of  time-dependent  loads. 


BACK-FACE  THERMAL  BOUNDARY  CONDITION 

As  stated  in  the  main  text,  the  analysis  of  the  temperature  history  given 
in  appendix  A  incorporates  the  assumption  that  no  heat  flow  takes  place  from 
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the  back  face  of  the  structural  layer .  This  assumption  follows  the  precedents 
set  by  previous  analyses  of  similar  configurations  and  leads  to  conservative 
calculated  temperatures.  The  purpose  of  this  section  is  to  provide  some  quan¬ 
tification  of  the  effects  of  that  assumption. 

In  order  to  carry  out  this  quantification,  a  finite-element  thermal  model 
was  formulated  by  using  the  SPAR  thermal  analyzer  (ref.  22).  In  the  SPAR  model, 
the  back  face  of  the  structural  layer  was  permitted  to  radiate  to  a  medium  at 
room  temperature.  Temperature  histories  were  calculated  for  the  eight  final 
designs  from  table  IV(d) .  Table  VII  contains  a  comparison  of  the  peak  tempera¬ 
tures  for  the  adiabatic  and  radiating  boundary  conditions.  The  largest  change 
in  peak  temperature  is  33  percent  for  titanium,  and  the  smallest  change  is  for 
Rene  41  and  carbon/carbon.  The  effect  of  radiation  might  have  been  expected 
to  be  largest  for  the  higher  temperature  materials.  Certainly  more  heat  is 
transferred  from  the  back  face  of  the  high-temperature  materials  than  the  low- 
temperature  materials.  However  the  difference  in  temperatures  from  analyses 
with  and  without  back-face  radiation  is  not  necessarily  greater  for  the  high- 
temperature  materials.  Other  considerations  are  involved  and  the  following 
explanation  is  offered.  Inspection  of  the  thermal  properties  and  the  response 
for  the  eight  materials  indicate  an  inverse  correlation  between  the  heat-storage 
dT2 

term  P2C2t2  -  and  the  difference  in  temperatures  in  the  two  columns  of 

table  VII.  Prcxn  elementary  energy  conservation,  a  high  value  of  heat  storage 
leads  to  a  high  rate  of  radiation  heat  transfer  from  the  back  face  of  the  struc¬ 
ture.  This  result  in  turn  requires  a  high,  structural  wall  temperature.  Con¬ 
versely,  the  materials  with  low-heat  storage  (notably  titanium  and  B/Al)  tend 
to  have  lower  wall  temperatures.  The  temperatures  in  the  left  column  of 
table  VII  are  upper  bounds  to  the  corresponding  temperatures  in  the  right 
column.  When  radiation  is  included,  the  higher  storage  materials  have  struc¬ 
tural  temperatures  closer  to  the  adiabatic  values  than  the  low-storage  mater¬ 
ials,  and  hence  the  effect  of  the  back  wall  radiation  is  greatest  for  these 
low-storage  materials. 

Although  the  effect  of  the  back-face  boundary  conditions  has  a  fairly 
sizable  effect  on  scane  peak  temperatures,  it  does  not  necessarily  have  a  large 
effect  on  design  mass.  As  an  example,  consider  the  largest  difference  in  tem¬ 
perature  being  that  of  titanium.  The  increase  in  allowable  stress  of  titanium 
from  548  K  to  451  K  (527°  F  to  352°  F)  is  about  16  percent.  If  the  structural 
mass  of  the  titanium  design  is  assumed  to  be  decreased  by  1 6  percent  and  the 
insulation  mass  remains  unchanged,  the  total  required  mass  is  decreased  by  about 
11  percent. 
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TABLE  I.-  THERMAL  PROPERTIES  USED  FOR  OPTIMIZATION  STUDIES  OF  INSULATED  PANELS 


TABLE  II,-  MECHANICAL  PROPERTIES  OF  STRUCTURAL  MATERIALS  USED  FOR  OPTIMIZATION  STUDIES  OF  INSULATED  PANELS 


(0.056)  (0.056)  (0.056)  (0.056)  (0.098)  (0.098)  (0.100)  (0.100) 


Titanium  Rene  41  Carbon/carbon  Lockalloy 


IftBLE  IV.-  CHARACTERISTICS  OF  MINIMUM-MASS  INSULATED  PANELS  FOR  VARIOUS  STRUCTURAL  MATERIALS 


(a)  Tpeak  =  1  089  K  (ISOOO  F) ;  N*  =  Ny  =  0;  N^y  =  0 


Material 

Total  mass 

Maximum 

temperature 

Time  of 
maximum 
temperature, 
sec 

Maximum 

stress 

ratio 

Time  of 
maximum 
stress  ratio, 
sec 

kg/m^ 

Ibm/ft^ 

K 

Op 

G  r  aph  i  te/po  ly  imi  de 

11.82 

2.42 

533 

500 

2000 

0 

— 

Graphite/epoxy 

15.87 

3.25 

450 

350 

2800 

.01 

0 

Boron/aluminum 

10.20 

2.09 

389 

600 

2000 

.10 

2000 

Aluminum 

16.06 

3.29 

450 

350 

2800 

0 

Titanium 

11.82 

2.42 

589 

600 

2000 

0 

Ren^  41 

6.30 

1.29 

1078 

1480 

800 

0 

Lockalloy 

7.82 

1  .60 

700 

800 

1600 

0 

Carbon/carbon 

1  .76 

.36 

1078 

1480 

800 

0 

(b)  Tpean  =  1089  K  (1500®  F);  N^  =  Ny  =  -525  kN/m  (-3000  Ibf/in.); 
Nxy  =  350  kN/m  (2000  Ibf/in.) 


Material 

Total  mass 

Maximum 

temperature 

Time  of 
maximum 
temperature, 
sec 

Maximum 

stress 

ratio 

Time  of 
maximum 
stress  ratio, 
sec 

kg/m^ 

Ibm/ft^ 

K 

Op 

G  r  aph i te/po ly imi de 

11  .96 

2.45 

533 

500 

2000 

0.65 

2000 

Gr  aph i t e/epoxy 

15.87 

3.25 

450 

350 

2800 

.87 

2800 

Boron/aluminum 

11  .13 

2.28 

583 

590 

1600 

1  .00 

1 600 

Aluminum 

17.58 

3.60 

450 

350 

2400 

1 .00 

2400 

Titanium 

13.18 

2.70 

589 

600 

2000 

.98 

2000 

Rene  41 

11.72 

2.40 

1078 

1480 

800 

1  .00 

800 

Lockalloy 

11.72 

2.40 

567 

560 

1600 

1  .00 

1600 

Carbon/carbon 

20.41 

4.18 

1081 

1486 

800 

1.00 

2400 

30 


I!^LE  rv.-  Continued 


(c)  Tpgai^  =  1089  K  (1500°  P) ;  Njj  =  Ny  =  -1.05  MN/m  (-6000  Ibf/in.); 
Nxy  =  700  kN/m  (4000  Ibf/in.) 


Material 

Total  mass 

Max imum 
temperature 

Time  of 
maximum 
temperature, 
sec 

Maximum 

stress 

ratio 

Time  of 
maximum 
stress  ratio, 
sec 

kg/m^ 

lbin/ft2 

K 

Op 

Graph i te/poly irai de 

13.28 

2.72 

533 

500 

2000 

0.76 

2000 

Gr  aph i t e/epoxy 

16.60 

3.40 

450 

350 

2800 

1.00 

2800 

Boron/aluminum 

14.84 

3.04 

522 

480 

1600 

1 .00 

1600 

Aluminum 

21.53 

4.41 

450 

350 

2000 

1  .00 

2000 

Titanium 

16.45 

3.37 

584 

592 

2000 

1.00 

2000 

Rene'^  41 

23.39 

4.79 

1078 

1480 

800 

1.00 

800 

Lockalloy 

16.55 

3.39 

503 

445 

1600 

1  .00 

1600 

Carbon/carbon 

41  .35 

8.47 

1076 

1477 

800 

1  .00 

2800 

(<3)  Tpeak  =  1089  K  (1500°  F) ;  N*  =  Ny  =  -1.58  MN/m  (-9000  Ibf/in.); 
Nxy  =  1.05  MN/m  (6000  Ibf/in.) 


Material 

Total  mass 

Maximum 

temperature 

Time  of 
maximum 
temperature, 
sec 

Maximum 

stress 

ratio 

Time  of 
maximum 
stress  ratio, 
sec 

kg/m^ 

lbm/ft2 

K 

Of 

G  r  aph  i  te/po  ly  imi  de 

13.96 

2.86 

530 

494 

n 

2000 

Gr  aphite/epoxy 

18.02 

3.69 

450 

350 

2400 

Boron/aluminum 

16.70 

3.42 

494 

429 

■RH 

2000 

Aluminum 

29.29 

6.00 

384 

232 

1.00 

2800 

Titanium 

20.36 

4.17 

548 

527 

1  .00 

2000 

Rene'^  41 

35.06 

7.18 

1078 

1480 

800 

1  .00 

800 

Lockalloy 

21  .53 

4.41 

475 

395 

1600 

1.00 

1600 

Carbon/carbon 

62.50 

12.80 

1018 

1373 

800 

1  .00 

2800 

31 


TABLE  IV. -  Continued 


(e)  Tpeak  =  1 089  K  (1500°  P) ;  =  Ny  =  -2.10  MN/m  (-12  000  Ibf/in.); 

^  Nvv  =  1.40  MN/m  (8000  Ibf/in.) 


Material 


Gr aph i te/poly imi de 

Gr aphite/epoxy 

Boron/aluminum 

Aluminum 

Titanium 

Rene"^  41 

Lockalloy 

Carbon/carbon 


Total  mass 

Maximum 

temperature 

kg/m^ 

Ibm/ft^ 

K 

Of 

17.53 

3.59 

498 

437 

20.07 

4.11 

440 

340 

19.48 

3.99 

499 

438 

33.88 

6.94 

384 

232 

24.17 

4.95 

570 

536 

46.72 

9.57 

1078 

1480 

26.80 

5.49 

456 

360 

82.61 

16.92 

924 

1204 

Time  of 
maximum 
temperature  i 
sec 


2000 

2400 

1600 

2400 

1600 

800 

1600 

800 


Maximum 

stress 

ratio 


1 .00 
.75 
1 .00 
1 .00 
1  .00 
1 .00 
1 .00 
1 .00 


Time  of 
maximum 
stress  ratio# 
sec 


2000 

2400 

1600 

2400 

1600 

800 

1600 

3600 


(f)  Tpeak  =  1089  K  (1500°  F) ;  =  Ny  =  -2.63  MN/m  (-15  000  Ibf/in.); 

Nxy  =  1.75  MN/m  (10  000  Ibf/in.) 


Material 

Total  mass 

Maximum 

temperature 

Time  of 
maximum 
temperature, 
sec 

Maximum 

stress 

ratio 

Time  of 
maximum 
stress  ratio, 

sec 

kg/ra^ 

Ibm/ft^ 

K 

Op 

G r aph i t e/po ly imi de 

21 .73 

4.45 

417 

290 

2000 

Gr  aph i t e/epoxy 

24.90 

5.10 

450 

350 

0 

Boron/aluminum 

22.36 

4.58 

517 

370 

1  .00 

2000 

Aluminum 

38.08 

7.80 

432 

318 

1.00 

2000 

Titanium 

28.12 

5.76 

555 

540 

1  .00 

1600 

Rene^  41 

58.44 

11.97 

1078 

1480 

1  .00 

800 

Lockalloy 

33.40 

6.84 

505 

449 

1600 

1  .00 

1600 

Carbon/carbon 

96.28 

19.72 

889 

1140 

1200 

1.00 

0 
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TABLE  rv.-  Continued 


(g)  Tpeak  =  1533  K  (2300°  F) ;  =  Ny  =  0;  N^y  =  0 


Material 

Total  mass 

Maximum 

temperature 

Time  of 
maximum 
temperature, 
sec 

Maximum 

stress 

ratio 

Time  of 
maximum 
stress  ratio, 
sec 

kg/m^ 

lbm/ft2 

K 

Op 

Gr aph i t e/po ly  im  ide 

15.77 

3.23 

533 

500 

2400 

0 

— 

Graphite/ epoxy 

19.53 

4.00 

450 

350 

3600 

.06 

0 

Bor on/a luminura 

13.57 

2.78 

589 

600 

2400 

.13 

2400 

Aluminum 

20.65 

4.23 

450 

350 

3600 

0 

— 

Titanium 

15.53 

3.18 

589 

600 

2400 

0 

— 

Ren^  41 

11.03 

2.26 

1089 

1500 

1200 

0 

— 

Lockalloy 

10.89 

2.23 

700 

800 

2000 

0 

— 

Carbon/carbon 

1  .76 

.36 

1515 

2268 

800 

.96 

0 

(•^)  Tpeak  =  1533  K  (2300°  F) ;  Njj  =  Ny  =  -525  kN/m  (-3000  Ibf/in.); 
Njjy  =  350  kN/m  (2000  Ibf/in.) 


Material 

Total  mass 

Maximum 

temperature 

Time  of 
maximum 
temperature, 
sec 

Maximum 

stress 

ratio 

Time  of 
maximum 
stress  ratio, 
sec 

kg/m2 

lbm/ft2 

K 

Op 

G r aph i t e/po ly im i de 

15.77 

533 

500 

2400 

0.26 

0 

Graphite/epoxy 

19.72 

450 

350 

3600 

.93 

3600 

Boron/aluminum 

14.89 

3.05 

560 

548 

2400 

1 .00 

2400 

Aluminum 

21  .43 

4.39 

450 

350 

3200 

1 .00 

3200 

Titanium 

16.40 

3.36 

589 

600 

2400 

1 .00 

2400 

Ren^  41 

15.28 

3.13 

1046 

1423 

1200 

1  .00 

1200 

Lockalloy 

13.38 

2.74 

656 

720 

1600 

1.00 

1600 

Carbon/carbon 

20.41 

4.18 

1520 

2277 

800 

1 .00 

0 
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TABLE  IV. “  Continued 


(i)  Tpeak  =  ’•533  K  (2300°  F) ;  =  Ny  =  -1,05  MN/m  (-6000  Ibf/in.); 

^  Nxy  =  700  kN/m  (4000  Ibf/in.) 


Material 


Graph ite/polyimide 

Graphite/epoxy 

Bor on/a lurainum 

Aluminum 

Titanium 

Rene*  41 

Lockalloy 

Carbon/carbon 


Total 

mass 

Maximum 

temperature 

kg/m^ 

Ibm/ft^ 

K 

_ 

o. 

16.21 

3.32 

500 

21  .04 

4.31 

332 

16.80 

3.44 

560 

549 

25.10 

5.14 

450 

350 

20.94 

4.29 

589 

600 

24.17 

4.95 

713 

824 

18.55 

3.80 

563 

554 

41  .35 

8.47 

1513 

2264 

Time  of 
maximum 
temperature 
sec 

2400 

4000 

2000 

2800 

2000 

2000 

2000 

800 


Maximum 

stress 

ratio 

Time  of 
maximum 
stress  ratio, 
sec 

0.94 

1 .00 

.91 

2000 

O 

o 

• 

2800 

.83 

2000 

1  .00 

2000 

1  .00 

2000 

1 .00 

0 

(j)  Tpeak  =  1533  K  (2300°  F) ;  Njj  =  Ny  =  -1.58  MN/m  (-9000  Ibf/in.); 
Njjy  =  1.05  MN/m  (6000  Ibf/in.) 


Material 

Total 

mass 

kg/m^ 

Ibm/ft^ 

Graph i te/poly imi de 

18.06 

3.70 

Gr  aphite/epoxy 

21  .58 

4.42 

Boron/aluminum 

19.24 

3.94 

Aluminum 

30.66 

6.28 

Titanium 

23.19 

4.75 

Rene*  41 

31  .64 

6.48 

Lockalloy 

23.63 

4.84 

Carbon/carbon 

62.40 

12.78 

Maximum 

temperature 


Time  of 
maximum 
temperature, 
sec 

Maximum 

stress 

ratio 

Time  of 
maximum 
stress  ratio, 
sec 

2000 

0.06 

0 

3200 

.95 

3200 

2000 

1  .00 

2000 

2800 

1  .00 

2800 

2000 

1  .00 

2000 

1600 

1  .00 

1600 

2000 

1  .00 

2000 

800 

1 .00 

0 

TABLE  IV. “  Concluded 


<k)  Tpeak  =  1533  K  (2300°  P) ;  =  Ny  =  -2.10  MN/m  (-12  000  Ibf/in.); 

Nxy  =  1.40  MN/m  (8000  Ibf/in.) 


Material 

Total  mass 

Maximum 

temperature 

Time  of 
maximum 
temperature, 
sec 

Maximum 

stress 

ratio 

Time  of 
maximum 
stress  ratio, 
sec 

kg/m2 

Ibm/ft^ 

K 

Op 

G  r  aph i te/poly imi de 

19.68 

4.03 

529 

492 

1 .00 

2000 

Gr aph it e/epoxy 

25-44 

5-21 

425 

1.00 

0 

Boron/aluminum 

22.02 

4.51 

499 

438 

1 .00 

2000 

Aluminum 

35-05 

7.18 

448 

346 

1  .00 

2400 

Titanium 

27.05 

5.54 

581 

586 

1  .00 

2000 

Rene^  41 

38.96 

7.98 

709 

817 

1600 

1.00 

1600 

Lockalloy 

28.51 

5.84 

474 

394 

2000 

1  .00 

2000 

Carbon/carbon 

82.56 

16.91 

1277 

1839 

1200 

1  .00 

0 

(1)  Tpeak  =  1533  K  (2300°  P) ;  =  Ny  =  -2.63  MN/m  (-15  000  Ibf/in.); 

Nxy  =  1.75  MN/m  (10  000  Ibf/in.) 


Material 

Total  mass 

Maximum 

temperature 

Time  of 
maximum 
temperature, 
sec 

Maximum 

stress 

ratio 

Time  of 
maximum 
stress  ratio, 
sec 

kg/m2 

lbm/ft2 

K 

Op 

Gr aph i te/poly imi de 

25.00 

5.12 

454 

358 

2400 

1 .00 

2000 

Gr  aphite/epoxy 

29.34 

6.01 

433 

319 

2800 

1.00 

0 

Boron/aluminum 

25.10 

5.14 

459 

367 

2400 

1.00 

2400 

Aluminum 

41  .06 

8.41 

412 

282 

2400 

1.00 

2400 

Titanium 

30.71 

6.29 

532 

498 

2000 

1.00 

2000 

Rene"*  41 

49.95 

10.23 

847 

1065 

1600 

1 .00 

Lockalloy 

33.79 

6.92 

463 

373 

2000 

1 .00 

2000 

Car  bon/carbon 

96.28 

19.72 

1220 

1736 

1200 

1 .00 

0 
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TABLE  V.-  OPTIMUM  MASS  OF  INSULATED  PANELS  FOR  VARIOUS  VALUES  OF 
HEATING  INTENSITY  AND  DURATION 

[Nx  =  Ny  =  -700  kN/m  (-4000  Ibf/in.);  Nxy  =  525  kN/m  (3000  Ibf/in.)] 


Material 

Tpeak  =  Toeg 

K  (1  500°  F) 

"^peak  ”  ^ 

K  (2300°  F) 

Tf  =  1500  sec 

Tf  =  2000  sec 

Tf  =  1500  sec 

Tf  =  2C 

100  sec 

kg/m^ 

Ibm/ft^ 

kg/m^ 

Ibm/ft^ 

kg/m^ 

Ibm/ft^ 

kg/m^ 

Ibm/ft^ 

Graph ite/polyimide 

11  .57 

2.37 

12.01 

2.46 

14.16 

2.90 

15.23 

3.12 

Graphite/epoxy 

15.14 

3.10 

15.97 

3.27 

18.31 

3.75 

20.16 

4.13 

Bor on/a lumi num 

*13.72 

*2.81 

15.14 

3.10 

*16.36 

*3.35 

19.33 

3.96 

Aluminum 

18.16 

3.72 

19.19 

3.93 

21  .43 

4.39 

23.00 

4.71 

Titanium 

*14.26 

*2.92 

14.45 

2.96 

16.80 

3.44 

18.46 

3.78 

Rene  41 

*16.70 

*3.42 

*16.75 

*3.43 

— 

Carbon/carbon 

*30.27 

*6.20 

*30.27 

*6.20 

*30.27 

*6.20 

*30.27 

*6.20 

Lockalloy 

*13.23 

*2.71 

*14.79 

*3.03 

*15.48 

*3.17 

*16.01 

*3.28 

*Teniperature  response  less  than  allowable  value. 
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TABLE  VI.-  EFFECT  OF  TIME-DEPENDENT  LOADS  ON  DESIGNS  OF  INSULATED  PANELS 


Z 

o^ 
II  00 
o 


>iW 

<J\ 
il  00 
o 


Xf-  Xr- 

23  ^ 

♦  M  *  II 

4c 


rO 


Eh 


nj 

Q) 

Pa 

EH 
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TABLE  VII.-  EFFECT  OF  BACK -FACE  RADIATION  ON  MAXIMUM  STRUCTURAL  TEMPERATURES 

Calculations  for  final  designs  corresponding  to: 

Nx  =  Ny  =  -1.58  MN/m  (-9000  Ibf/in.); 

Nxy  =  1.05  MN/ra  (6000  Ibf/in.);  Tf  =  2000  sec; 

Tpeak  =  1089  K  (1500°  F) 


Material 

Maximum  temperature 

Adiabatic 

back  face 

Radiating 

back  face 

K 

Op 

K 

Op 

Graph ite/polyimide 

530 

494 

453 

355 

Graphite/epoxy 

450 

350 

393 

248 

Boron/aluminum 

492 

429 

441 

334 

Aluminum 

384 

232 

362 

191 

Titanium 

548 

527 

451 

352 

Rene  41 

1078 

1480 

1069 

1464 

Lockalloy 

475 

395 

448 

346 

Car bon/ car bon 

1018 

1373 

963 

1274 
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Operating  temperature 

Figure  3.-  Relationship  between  operating  temperature  and  mass  of 

insulated  panel. 


Temperature 


Time 

Figure  4.-  Discrete-times  approach  for  constraints  on  transient 

temperature  response. 


Number  of  discrete  times 


Figure  5.-  Effect  of  number  of  discrete  times  on  optimum  ma^s  of 
insulated  Gr/E  panel.  =  Ny  =  -700  kN/m  (-4000  Ibf/in.); 
Nxy  =  525  kN/m  (3000  Ibf/in.);  Tpeak  =  K  (1 500^  f) ; 

Tf  =  1500  sec. 


Temperature  or  stress  ratio 


•e  or  stress  ratio 


Number  of  discrete  times 


Figure  7 Effect  of  number  of  discrete  times  on  optimum  mass  of 
insulated  Gr/E  panel.  Nx  =  Ny  =  -1.05  MN/m  (-6000  Ibf/in.); 
Nxy  =  700  kN/m  (4000  Ibf/in.);  Tpggi^  =  1089  K  (1500°  F) ; 

=  1500  sec. 


Optimum  mass,  kg/m 


Carbon/carbon 


Rene  41 


Aluminum 


Lockalloy 


^  Titanium 
'^Gr/E 

^Xb/ai 


Rene  41 
Lockalloy 


B/Al 


Carbon/carbon 


-15  X  10' 


Ibf/in 


Figure  8 


0  :i:0  -1.5  -2.0  -2.5 

Load,  MN/m 

Optimum  mass  of  insulated  panels.  Ny  =  Njj/  Nxy  “  (2/3) Njj; 
=  1089  K  (1500°  F)  . 


Carbon/carbon 


0  -.5  -1.0  -1.5  -2.0  -2.5  -3.0 

Load,  MN/m 

-  Optimum  mass  of  insulated  panels.  Nx  =  Ny;  Nxy  =  (2/3) N^; 

=  1533  K  (2300°  F)  . 


Optimum  mass,  kg/m 


Ibm/ft 


Edge -fixity  coefficient,  )3 


Figure  11.-  Effect  of  restrained  thermal  expansion  on  mass  of 
insulated  panels.  =  Ny  =  -1.58  MN/m  (-9000  Ibf/in.); 
Nxy  =  1.05  MN/m  (6000  Ibf/in.);  Tpeak  =  K  (1500°  F)  ; 

T4=  =  2000  sec. 


Load  factor,  F  (t) 
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